Quantum gravity effects on compact star cores 
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Abstract 



Using the Tolman-Oppenheimer-Volkoff equation and the equation of state of zero temperature 
ultra-relativistic Fermi gas based on GUP, the quantum gravitational effects on the cores of compact 
stars are discussed. Our results show that 2m{r)/f varies with f. Quantum gravity plays an important 
Q^l role in the region f ~ 10'^. Furthermore, near the center of compact stars, we find that the metric 
^ ' components are gu ~ r'^ and g^r = (1 — All these results are different from those obtained 

from classical gravity. 
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\ The configuration of a spherically symmetric static star, composed of perfect fluids, is determined 

by the Tolman-Oppenheimer-Volkoff (TOV) equation (in c.g.s. units) [H [2] 



dr yi"^ l"^ > r[r-2Gm{r)/c^] ' ^ ' 

^ ' with 

S : ^ = 4vr.V(r), (2) 

dr 

where c is the velocity of light. G is the gravitational constant. P and p are respectively the pressure 
and the macroscopic energy density measured in proper coordinates. Supplied with an equation of 
state and appropriate boundary conditions, eqn. ([T]) and eqn. ([2]) determine P{r), m{r) and p{r). If 
the pressure and gravitational potential is everywhere small, i.e., P{r) <^ pc^, 2Gm{r)/c^r <C 1, the 
TOV equation reduces to the fundamental equation of Newtonian astrophysics 
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Most of the low density compact stars like white dwarfs are well described by Newtonian gravity. 
However, for compact stars like neutron stars and other exotic compact stars, general relativity plays 
an important role [3]. Although the real structures of neutron stars and other Fermi stars are very 
complicated, the most important physical property is given by a simple model in which the pressure 
of cold degenerate fermions contends against the gravitational collapse [2]. 

For the high density and high pressure cold Fermi stars it is worth considering the effects of 
quantum gravity. In the absence of a full theory of quantum gravity, effective models are useful 
tools to gain some features from quantum theory of gravity. One of the most important models is 
the generalized uncertainty principle (GUP), derived from the modified fundamental commutation 
relation H El EJ [3 [HJ El [TO] 

[x,p\=ih{l + pp'), (4) 

where /3 = Po§/h'^ = Po/c^M^, = Gh/c^, = hc/G. h = h/lir is the Planck constant and /3o 
is a dimensionless parameter. With this modified commutator, one can easily derive the generalized 
uncertainty principle (GUP) 

AxAp>^[l + /3{Ap)% (5) 

which in turn gives the absolutely smallest uncertainty in positions, i.e., the minimum measurable 
length 

Ax > Ainin = = V^lp- (6) 

Note that the model in ^ considers only the minimal uncertainty in position. In this case, the 
quantum mechanics structure underlying the GUP has been studied in full detail [7]. The statistics 
of ideal gases based on GUP has been discussed by many authors [111 [T2l \T3\ [T4l I15j . In our recent 
work, we have studied a system composed of zero temperature ultra-relativistic Fermi gas based on 
GUP [15]. The Newtonian equation with uniform pressure was employed to discuss stellar structures. 
In this paper, we will apply the TOV equation to investigate the cores of the stars where quantum 
gravity plays a leading role. 

In [15], the particle number, energy density and pressure for an ultra-relativistic system were given 

V - (^-^/(^)' 

" - -^E%hi.), (8) 
Stt 



c^{hc 

EUi^), (9) 



(he) 

where Eh = c/^/P = MpC^/y^ denotes the Hagedorn energy, introduced in \X5l and k = ep 



2 



Bp /Eh- Moreover 



9{k) = 



4(1 + ^2)2' 



1 



k{k? - 1) _i , . 



(1 + K^) 
Kf{K) - h{K). 



(10) 
(11) 

(12) 



It is worth noting that when k increases, the pressure blows up, while the energy density and number 
density are both bounded. This is a manifestation of the minimal length. 

Two configurations of compact stars have been addressed in [15], by applying the Newtonian limit 
eqn. (j3]) with uniform density. One is that the star is almost composed of ultra-relativistic particles. 
The other is that the major contribution to the mass is from non-relativistic cold nuclei. However, 
to discuss the core of ultra-compact stars like neutron stars, one should use TOV equations Setting 
r = rgf, m = mom, P = PqP and 

Gmo 



Pop, Pq = Poc 



1, 



the TOV eqn.s ([T|) and ^ are reduced to the following dimensionless ones 

dP _ ^~ j_ 'fn + P 
df r{f — 2m) ' 



dm 
df 



f'^p. 



(13) 
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When there is no introduction of quantum gravity, for a system almost composed of ultra-relativistic 
fermions, the equation of state is P = p/3. An exact solution is given in |16j 



2m(f) 3 ~ 1 

= ^, P{r) = —r \ 
r 7 14 



(16) 



The pressure is not zero on the surface of the star. This does not meet the physical boundary 
conditions. However, the point is that it is an analytic solution describing the central region of 
compact stars with divergent pressure in the center [2j. Note that the length scale ro in eqn. (I13p is 
uncertain. Thus r, m, p and P can be any size. 

From eqn. (jl6p . the pressure is divergent in the center. Therefore, influences from quantum gravity 
should be included in the discussion. Obviously, near the surface, particles are non-relativistic while 
in the region around the center, particles are ultra-relativistic |2j. This determines the equations of 
state and boundary conditions. In the vicinity of r = 0, the equation of state is given by eqn. ([7]), 
eqn. ([8]) and eqn. ([9]). Under the limit k — )• 0, it is straightforward to recover P = p/Sc^. Defining 
r = rof, m = mom with 
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(17) 
(18) 
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where ro is the minimum radius in [15] 



r-o = y ^/3o/p - V^A^in. (20) 

The expressions (jl9p and (j20p show the system can not be arbitrary scale, determined entirely by 
This indicates that our discussion is focused on the central region of compact stars. Substituting the 
above expressions for P and p (eqn. ([8]) and eqn. ([9])) into eqn. (jl4p and eqn. (jl5p . one gets 



drh{rj ,2 



(iK(f) — K(f) [m(f) + f^g{K)\ 



f^h{K), (21) 

(22) 



df f [f — 2m(f )] 

Since the density is regular in the center, one has m(0) = as a boundary condition. After setting 
kq = k(0) as another boundary condition, eqn. ()2ip and eqn. ()22p are integrated numerically in Table 
I - Table V. 

In Table I to Table III, we perform the integration with different K{r). Four conclusions can be 
drawn from these tables: 

• Different from the results obtained in classical gravity, 2m(f)/f varies with f but not a constant 
3/7. For example, with K(f) = 0.1 in Table I, the deviation of 2m(f)/f is about 4%. 

• 2m{r)/f is not sensitive to different initial value kq. 

• For large K(f) or small f, quantum gravity contribution is important to the value of 2m{f)/f. 
As K(f) decreases, or f increases, the configuration approaches the classical one obtained in [16j . 
with a constant 2m{f)/f = 3/7. 

• Quantum gravity plays an important role in the region r ~ lO^ro. 
Some analytic solutions can be obtained in extreme cases as follow. 

• Under k — )• 0, it is easy to see that /i(k) ~ k^/4, g{K) ~ Then from eqn. (I2ip and eqn. 
(j22]) . we obtain 

^"^'^ - ^, <r) = (^)'^^-^/^ P{r) = = 3^f-^ for large f. (23) 



1/4 

r 

This solution is nothing but the classical one without quantum gravity. 



• Under r — )• and k — )• 00, eqn. (j2T|) and eqn. (p2|) can by replaced by asymptotic expressions 

(24) 



dm{f) 1 _2 



df 4 



4 



dr 



-K[r 



r\r 



2m{r)] 



The solution of these equations is 



rh{r) 



12 



, K{f) 



32 1 



, P(f) 



for f — )• 0. 



(25) 



(26) 



The solution (j26p represents the situation where quantum gravity dominates. This happens near the 
center of ultra-compact stars. One can see that it is quite different from the solution of classical 
gravity. Table IV is the numerical result integrated for large ^(f), well consistent with the asymptotic 
solution (j26]) . 

For a spherically symmetric static compact star, the metric is given by [3] 



Mr) 



1 



2Gm{r) 



9tt 



, , 1 dB 
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m(r) + 





2m{r 
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r 






1 - 


2Gm 


C2 




c^r 



-1 



Then for r — >• 0, from (1261). we have 



(27) 



(28) 



A(r) 



1 -f2/6' 

One may compare ()29p with the classical results 



B{r) ~ P. 



A(r) 
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4' 



B{r) 



(29) 



(30) 



In Table V, eqn. (j2ip and eqn. (|22l) are integrated with a large initial kq- It is interesting that 
2m/f reaches a maximum value 0.734 in the vicinity of r = 3.00 rg. Our calculation shows that near 
the center, 2m/f = f2/6 which indicates that 2rh/r increases with r. On the other hand, as f — t- c«, 
2rh/r — )• 3/7. Therefore, the maximum of 2rh/f at r = 3.00 tq is a turning point, where quantum 
gravity effect starts to dwindle. From Table V, one also finds that has a small range of fluctuation. 
A minimum (1 - 0.279)"^ = 1.39 is achieved at r ~ 12.5 tq. This minimum is about one-third of the 
maximum (1 — 0.734)^^ = 3.76 at r ~ 3.00 tq. We do not have good explanation for this fluctuation. 
It may be caused by the effectiveness of our model. Finally, tends to the constant 7/4 at large r 
as expected. The profile of 2m/f versus r is plotted in Fig[TJ One can see that the upper limit, 8/9 
on the surface of a spherically symmetric static star, is well satisfied. 

In summery, we discussed the structure of ultra-compact star cores by a simple effective quantum 
gravity model. The model, GUP, introduces a new equation of state, determined by eqn.s d?]), ([8]) and 
([9]). By plugging the equation of state into TOV equations, we found some different features from 
previous works in literature. 

Since quantum gravitational effects play an important role only in high density, we considered 
configurations in which a star is almost composed of ultra relativistic particles. The asymptotic 
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Kq 




m(r) 


f 


2m(r)/r 

\ /I 


1000.0 




21.88 


97.70 


0.448 


100.0 




21.87 


97.70 


0.448 


50.0 




21.85 


97.70 


0.447 


20.0 




21.79 


97.60 


0.446 


10.0 


0.1 


21.71 


97.60 


0.445 


8.0 




21.65 


97.50 


0.444 


5.0 




21.53 


97.50 


0.442 


3.0 




21.29 


97.40 


0.437 


1.0 




19.72 


95.20 


0.414 


0.5 




18.36 


87.90 


0.418 



Table I: Integration from kq to K,(f) = 0.1. The value of 1m,(f)lf is insensitive to initial condition kq- 1fn(f)lr 
has relatively large deviation from 0.429 since at f = 97.7 quantum gravity has evident effects. 





K(f) 


m(f) 


f 


2rh{r)/r 


1000.0 




1985.47 


9250.70 


0.429 


100.0 




1985.48 


9250.80 


0.429 


50.0 




1985.51 


9251.00 


0.429 


20.0 




1985.61 


9251.60 


0.429 


10.0 




1985.76 


9252.50 


0.429 


8.0 


0.01 


1985.86 


9253.00 


0.429 


5.0 




1986.14 


9254.40 


0.429 


3.0 




1986.72 


9257.20 


0.429 


1.0 




1988.16 


9269.70 


0.429 


0.5 




1981.62 


9265.40 


0.428 


0.1 




2016.18 


9405.60 


0.429 



Table II: Integration from kq to ^(f) = 0.01. The different initial value kq has almost no effect on 2m(f)/f. f 
is large enough to overwhelm quantum gravity influences. 

solutions near the center are given by (j26p and (|29p . The complete picture is given by numerical 
calculation. Quantum gravitational effects play a leading role only in a relatively small range ~ 
103 ro = 10^ VA) ^min- Outside this region, the solutions are determined by eqn (jl6p and ()30p . Our 
discussion can be applied to neutron stars, for example. 
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Kq 




m(f) 


r 


2rn,(f) /r 


1 r\r\r\ r\ 

lUUU.U 




198o72.71 


nocTT^o err* 

925778.50 


rt A on 

0.429 


iUU.U 




190072.00 


925778.70 


n /I on 


olJ.L) 




198373.07 


925778.90 


n /I on 

0.429 


on n 
zU.U 




iyOo ( o.Do 


y2o ( ( y.oU 


n /I on 


10.0 


0.001 


198374.58 


925780.30 


0.429 


5.0 




198376.38 


925781.40 


0.429 


1.0 




198393.76 


925802.80 


0.429 


0.1 




198516.29 


925868.20 


0.429 


0.01 




201594.58 


940268.10 


0.429 



Table III: Integration from kq to K{f) — 0.001. The numerical results completely match the asymptotic solution 
eqn. 





K;(f) 


m{f) 


f 


2rh[f)/f 


1000.0 


20 


2.85 X 10-^ 


0.700 


8.15 X 10'^ 


500.0 


2.77 X 10"^ 


0.693 


7.99 X 10"^ 


200.0 


2.52 X 10"^ 


0.672 


7.51 X 10"^ 


100.0 


2.14 X 10"^ 


0.636 


6.72 X 10"^ 


50.0 


1.41 X 10"^ 


0.554 


5.10 X 10"^ 


30.0 


0.60 X 10"^ 


0.417 


2.89 X 10"^ 



Table IV: Integration from kq to K{f) — 20. The large value of n corresponds to r — > 0. 2rfi{f)/f depends 
sensitively on f. 



Kq 


K(f) 


rh{f) 


f 


2rh{r)/f 




100.0 


2.32 X 10"^ 


0.303 


1.53 X 10"^ 




50.0 


7.05 X 10-^ 


0.439 


3.21 X 10"^ 




20.0 


2.85 X 10"^ 


0.700 


8.14 X 10"^ 




10.0 


7.87 X 10"^ 


0.984 


0.160 




5.0 


0.205 


1.365 


0.301 


1000.0 


1.0 


0.983 


2.706 


0.727 




0.8 


1.100 


3.001 


0.734 




0.7 


1.169 


3.220 


0.726 




0.5 


1.344 


4.052 


0.663 




0.3 


1.825 


8.091 


0.451 




0.25 


2.359 


12.465 


0.279 




0.2 


4.062 


22.801 


0.356 




0.1 


21.880 


97.689 


0.448 




0.01 


1985.47 


9250.70 


0.429 



Table V: Integration with a fixed kq = 1000. 2m/f reaches its maximum value 0.734 in the vicinity of 
r = 3.00 tq. Irnjr has a small range of fluctuation and achieves a minimum value 0.279 at r ~ 12.5 ro. 
Eventually, Ifnjf tends to the constant 7/4 at large f. 



7 



•••• 



Figure 1: For a fixed kq = 1000, 2m/r versus the radius r = f. As f — t- 0, 2m/r ~ r^. 2m/r has 
a maximum around f = 3. 2m/r acquires the asymptotic value 0.429 at large r. The dashed line 
represents 2m/ r = 0.429 while the dotted line represents 2m/ r = 8/9, the upper limit of 2m/ r on the 
surface of a spherically symmetric static star. 
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